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Defects  created  during  the  manufacture  of  an  oxide/oxide  and  two  non-oxide  (SiC/SiNC  and  MI  SiC/SiC) 
ceramic  matrix  composites  (CMCs)  were  categorized  as  follows:  (1)  Intra-yarn  defects  such  as  dry  fibers, 
(2)  Inter-yarn  defects  such  as  those  at  crossover  points,  matrix  voids,  shrinkage  cracks  and  interlaminar 
separation,  and  (3)  Architectural  defects  such  as  layer  misalignment.  Their  impact  on  elastic  properties 
was  analytically  investigated  using  a  stiffness  averaging  approach  considering  the  defects  to  have  volu¬ 
metric  and  directional  influences.  In-plane  tensile  and  shear  moduli  as  well  as  the  through-thickness 
compressive  modulus  were  experimentally  evaluated.  Results  of  analytical  model  were  around  7%  on 
average  from  the  mean  value  of  the  experimental  data.  It  was  observed  that  interlaminar  separation  dras¬ 
tically  reduced  the  through-thickness  modulus  by  about  63%  for  the  SiC/SiNC,  40%  for  the  MI  SiC/SiC  and 
around  32%  for  the  oxide/oxide  composites.  Shrinkage  cracks  in  oxide/oxide  composite  reduced  the  in¬ 
plane  tensile  and  shear  moduli  by  14%  and  8.8%,  respectively. 

©  2013  Elsevier  Ltd.  All  rights  reserved. 


1.  Introduction 

A  multitude  of  experimental  investigations  have  revealed  the 
vulnerability  of  CMCs  to  strength  degradation  and  damage  upon 
prolonged  exposure  to  service  loads  and  environmental  conditions, 
especially  at  a  high  level  of  porosity  and  in  the  presence  of  cracks 
typically  created  during  manufacture  (see  for  example,  [15]).  Cur¬ 
rent  design  and  modeling  approaches  have  been  unable  to  improve 
upon  our  understanding  of  the  mechanisms  responsible  for  the  on¬ 
set  of  life-limiting  damage  in  the  presence  of  defects.  Such  lack  of 
insight  has  thus  far  limited  the  acceptance  of  CMCs  in  industrial 
applications  and  without  the  ability  to  identify  defects  and  quan¬ 
tify  their  effect  on  the  properties  of  CMCs,  further  development 
of  CMC  liners  and  other  engine  parts  would  have  to  rely  largely 
on  expensive  trial  and  error  approaches  involving  high  tempera¬ 
ture,  long  term  testing.  Besides  being  expensive,  the  trial  and  error 
approach  may  produce  components  with  uncertain  structural  life 
and  may  not  fully  exploit  all  the  desirable  characteristics  of  CMC. 
Thus  a  greater  fundamental  understanding  of  the  impact  of  defects 
on  the  behavior  of  CMCs  is  of  crucial  importance  for  further 
development  of  these  materials. 
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A  few  research  activities  reported  in  archived  literature  were 
devoted  to  understand  the  distribution  of  defects  and  their  impact 
on  the  thermal  and  mechanical  properties  of  CMC.  Del  Puglia  et  al. 
developed  a  strategy  to  model  the  effect  of  defects  on  thermal 
transport  properties  of  a  CMC  [4,6].  Physical  defects,  such  as  matrix 
voids  and  cracks,  were  identified  and  categorized  and  the  effect  of 
each  category  on  the  thermal  transport  through  the  composite  was 
numerically  modeled.  The  composite  under  investigation  was  a  10 
layer  T300  carbon/SiC  composite  in  which  carbon  fabric  was 
impregnated  using  a  polymer,  burned  in  air,  then  infiltrated  with 
silicon.  The  silicon  reacted  with  the  carbon  residues  forming  silicon 
carbide.  The  authors  proposed  four  types  of  defects  for  this  com¬ 
posite  system:  (A)  inter-fiber  micro-porosity  including  spherical 
voids  and  cracks  parallel  to  fiber  tows,  (B)  trans-tow  cracks  run¬ 
ning  through  the  tows  in  the  planes  parallel  to  the  fibers,  (C)  ma¬ 
trix  cracks  perpendicular  to  the  direction  of  the  tows  and  (D) 
matrix  voids  which  occur  at  the  intersection  of  orthogonal  tows. 
Del  Puglia  et  al.  concluded  that  the  dominant  defects  types  are  A 
and  B  and  that  both  considerably  degrade  the  thermal  properties 
of  the  composite  [5]. 

Computed  tomography  (CT)  scanning  was  used  to  characterize 
matrix  voids  in  a  Ml  SiC/SiC  composite  and  finite  element  analysis 
was  implemented  to  calculate  the  localized  stress  field  [1]. 
Multiple  images  of  slices  from  a  single  specimen,  with  an  average 
void  content  of  6-7%,  were  used  to  render  a  3D  image  containing 
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Nomenclature 

CMC 

ceramic  matrix  composites 

rf 

fiber  radius 

SiNC 

silicon-nitride-carbide 

Sf 

fiber  strength 

E 

elastic  modulus 

T 

interfacial  shear  strength  between  the  fiber  and  the  ma¬ 

G 

shear  modulus 

trix 

K 

bulk  modulus 

L 

distance  between  defects 

V 

Poisson’s  ratio 

Av 

vector  of  area  fraction  of  defects 

Subscripts 

P 

distance  between  yarns 

fl 

fiber  in  longitudinal  direction 

a,  b 

long  and  short  diameters  of  the  yarn  elliptical  cross-sec¬ 

d 

damaged 

tion 

w 

warp  yarns 

dz 

defect  dimension  in  z  direction 

f 

fill/weft  yarns 

df 

defect  dimension  in  fill/weft  direction 

m 

matrix 

dw 

defect  dimension  in  warp  direction 

V 

voids 

Vf 

fiber  volume  fraction 

porosities  and  defects.  FEA  using  ANSYS®  and  ABAQUS®  were  uti¬ 
lized  to  create  a  mesh  and  evaluate  the  effect  of  the  defects  while 
modeling  the  composite  as  an  orthotropic  homogenous  material. 
The  study  concluded  that  small  voids  had  a  little  impact  on  Young’s 
modulus  but  caused  an  increase  in  localized  stress  and  thus  had  an 
effect  on  the  onset  of  cracks. 

Peters  et  al.  investigated  the  effect  of  voids  on  the  elastic  moduli 
of  SiC/SiC  composites  made  from  Nicalon  fibers  coated  with  a  CVl 
pyro-carbon  coating  and  a  CVl  matrix  [20].  The  void  volume  frac¬ 
tion  was  estimated  to  be  around  14  ±  3%.  Peters  et  al.  utilized  the 
model  developed  by  Hashin  [12]  to  calculate  the  effect  of  voids 
on  the  longitudinal  and  shear  moduli  as  well  on  the  transverse 
moduli  of  the  composite.  The  effect  of  voids  was  included  in  the 
calculation  of  the  composite  properties  by  either  reducing  the 
modulus  of  the  matrix,  or  reducing  its  volume  fraction.  They  con¬ 
cluded  that  there  is  little  influence  of  the  voids  on  the  in-plane  lon¬ 
gitudinal  tensile  modulus  while  a  major  impact  was  witnessed  on 
the  shear  and  transverse  in-plane  moduli  of  up  to  60%  as  compared 
to  the  theoretical  values  of  fully  dense  materials. 

Liaw  et  al.  analytically  studied  the  effect  of  voids  for  metal  and 
ceramic  matrix  composites  by  considering  the  composite  to  be  an 
infinite  solid  plate  with  periodically  distributed  inhomogeneities 
[14].  The  homogeneous  average  strains  along  the  solid  have  peri¬ 
odic  disturbances  due  to  the  existence  of  defects.  Liaw  et  al.  used 
Eshelby’s  theory  to  calculate  the  strain  distribution  in  the 
composite. 

It  can  be  seen  that  two  modeling  approaches  emerge  out  of  the 
literature  review  presented  above: 

(1 )  A  homogenization  approach  in  which  the  composite  is  mod¬ 
eled  as  a  homogeneous  anisotropic  body  with  periodic  [14] 
or  non-periodic  [1,20]  disturbances. 

(2)  A  meso-volume  based  approach  in  which  the  composite  is 
divided  into  micro-subcells.  Each  subcell  includes  a  part  of 
the  preform  with  or  without  defects  [5]. 

The  homogenization  approaches  viewed  the  composite  as  a 
homogenous  anisotropic  body  and  modeled  the  impact  of  defects 
as  disturbances  to  its  homogeneity,  neglecting  the  effect  of  the 
location  and  size  of  defects,  with  respect  to  the  fabric  preform. 
On  the  other  hand,  the  meso-volume  approach  accounted  for  the 
location  and  size  of  defects  as  they  relate  to  the  yarn  for  a  specific 
type  of  CMC  and  modeled  their  impact  on  the  thermal  transport 
properties  using  EEA.  It  correlated  defects  to  a  single  yarn,  but 
did  not  provide  an  understanding  of  the  defect  distribution  and 
shape  with  respect  to  the  fabric  preform. 


In  this  work,  defects  are  categorized  based  on  their  location 
with  respect  to  the  fabric  preform  and  the  yarns  for  three  CMCs 
and  used  to  model  the  effect  of  defects  on  their  elastic  properties. 

2.  Materials 

Three  composite  systems  were  investigated  in  this  work: 

(1)  SiC/SiNC  composite  comprising  a  non-stochiometric  SiC  (CG 
Nicalon™)  fiber,  coated  with  boron  nitride,  in  a  matrix  of  sil¬ 
icon,  nitrogen  and  carbon  manufactured  by  multiple  itera¬ 
tions  of  a  polymer  pyrolysis  process.  The  fiber  preform  was 
a  cross-ply  balanced  8  harness  satin  weave  with  24  ends 
per  inch  and  42%  fiber  volume  fraction. 

(2)  Oxide/Oxide  composite  comprising  a  Nextel™  720  AI2O3 
fiber,  without  coating,  in  a  sol  gel  matrix  of  Alumina  Silicate. 
This  composite  has  a  nano-porosity  of  25%  in  the  matrix  for 
crack  deflection.  The  fiber  preform  was  a  cross-ply  balanced 
8  harness  satin  weave  with  27  ends  per  inch  and  46%  fiber 
volume  fraction. 

(3)  Melt  Infiltrated  in  situ  BN  SiC/SiC  composite  comprising  a 
stochiometric  SiC  (Sylramic™)  fiber,  with  an  in  situ  boron 
nitride  treatment  (performed  at  NASA-GRC),  and  coated 
with  a  thin  layer  of  BN.  The  SiC  matrix  is  infiltrated  by  vapor 
deposition  followed  by  slurry  casting  of  SiC  particulates  and 
a  final  melt  infiltration  of  Si.  The  fiber  preform  was  a  cross- 
ply  balanced  5  harness  satin  weave  with  20  ends  per  inch 
and  36%  fiber  volume  fraction. 

3.  Experimental  evaluation  of  mechanical  properties 

Three  different  elastic  properties  were  experimentally  evalu¬ 
ated  for  the  composite  systems  listed  above: 

(1)  In-plane  elastic  tensile  modulus  was  measured  using  tensile 
tests  performed  per  NASA’s  Enabling  Propulsion  Materials 
(EPM)  testing  standards  (equivalent  to  ASTM  Cl 359).  The 
elastic  modulus  was  calculated  as  the  slope  of  the  stress- 
strain  curve  in  the  linear  region.  Ml  SiC/SiC  and  oxide/oxide 
were  tested  at  room  temperature  and  SiC/SiNC  was  evaluated 
at  temperatures  ranging  from  room  temperature  to  1100  °C. 

(2)  In-plane  shear  modulus  was  determined  using  biaxial 
extensometry  on  samples  machined  out  of  the  panels  at 
45°  and  shear  analysis  consistent  with  ASTM  D3518  at  room 
temperature.  The  curves  were  fitted  in  the  linear  range  of 
the  stress  strain  curve  to  determine  the  modulus. 
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(3)  Through-thickness  compressive  elastic  modulus  was  evalu¬ 
ated  at  room  temperature  using  a  series  of  stacked  disks  in 
conformance  with  the  technique  developed  by  Ojard  et  al. 
[17].  In  this  experiment,  each  individual  disk  was  ground  flat 
to  remove  asperities  and  enough  disks  were  machined  for  a 
2.54  cm  extensometer  to  be  flagged  onto  the  sample. 
Through  the  center  of  each  disk  a  hole  was  machined  so  that 
a  graphite  rod  could  be  inserted  to  hold  the  stack  in  place  and 
eliminate  disk  movement  during  initial  loading.  The  rod  was 
machined  short  so  that  it  would  not  sustain  any  load  that 
could  affect  the  measured  strain.  Even  though  the  disks  were 
machined  flat,  there  is  a  typical  initial  compliance  to  the 
stack  that  had  to  be  overcome  by  sufficient  load.  Once  the  lin¬ 
ear  region  of  the  stress-strain  curve  was  reached,  stress 
strain  data  were  fitted  to  determine  the  modulus. 

4.  Geometric  categorization  of  defects 

Micrographic  images  of  pristine  samples  were  taken  for  the 
three  composite  systems  under  consideration  from  the  same  plates 
that  were  mechanically  tested.  Areas  without  a  matrix  material  or 
with  distorted  fabric  architecture  were  considered  as  defects.  It  is 
important  to  note  that  each  one  of  these  images  only  provide  a  pla¬ 
nar  view  of  the  three  dimensional  defects.  Based  on  these  micro¬ 
graphic  images,  defect  types  were  categorized  as  follows: 

(1)  Intra-yarn  defects  that  are  located  inside  the  yarn  and  lim¬ 
ited  by  the  outer  layer  of  fibers  defining  the  yarn.  Only  one 
type  of  inter-yarn  defect  was  identified: 

a.  Dry  fibers:  In  this  defect  a  percent  of  the  fibers  in  the  yarn 
are  not  coated  with  matrix.  This  defect  can  occur  due  to 
reasons  related  the  yarn  itself  such  as  in  the  case  of  a 
yarn  with  a  high  packing  factor.  Such  defect  can  also 
occur  due  to  reasons  related  to  the  fabric  preform  when 
matrix  infiltration  path  is  obstructed  by  another  yarn 
within  the  same  layer  or  from  a  different  layer.  Fig.  1 
shows  a  micrographic  image  of  a  SiC/SiNC  composite 
along  with  a  planar  view  showing  a  possible  location 
for  this  type  of  defect. 

(2)  Inter-yarn  defects  that  are  not  limited  to  the  yarn  volume 
and  may  occur  either  outside,  or  inside  and  outside  the  yarn 
at  the  same  time.  Four  different  subcategories  were 
identified: 

a.  Defects  at  crossover  points:  These  defects  are  in  the  form 
of  voids  located  under  the  crimped  portions  of  the  yarn 
at  crossover  points.  Their  width  is  limited  by  the  width 
of  the  yarn  and  their  height  by  the  orthogonal  yarn  at 


Fig.  1.  Dry  fibers  in  a  micrographic  image  and  as  represented  using  pcGINA©.  The 
dark  object  located  inside  the  yarn  in  the  figure  on  the  right  represents  this  type  of 
defect. 


the  crossover  point.  These  defects  are  highly  affected 
by  the  architecture  of  the  preform.  For  example,  there 
are  only  25  crossover  points  in  a  unit  cell  of  a  5  harness 
satin  weave,  of  which  only  five  locations  have  the  possi¬ 
bility  of  forming  defects  because  this  is  where  yarns 
crimp  to  cross  over  other  yarns.  There  are  four  sides  to 
each  crossover  location,  so  the  maximum  number  of  pos¬ 
sible  defects  for  this  type  is  5  x  4  =  20.  These  defects  are 
caused  due  to  the  obstruction  of  the  crimped  yarn  to  the 
flow  of  the  matrix  material.  Fig.  2  shows  a  micrographic 
image  of  a  SiC/SiNC  composite  along  with  planar  and  side 
views  showing  possible  locations  of  this  type  of  defects. 

b.  Matrix  voids:  Pockets  of  air  and/or  other  gases  can  be 
entrapped  in  the  matrix  during  infiltration  causing  voids. 
These  voids  tend  to  be  spherical  in  shape  and  to  exist  in 
the  matrix  material  outside  the  yarns.  Fig.  3  shows  pos¬ 
sible  locations  of  this  type  of  defects  in  a  micrographic 
image  and  a  side  view. 

c.  Shrinkage  cracks:  This  type  of  defect  occurs  when  the 
matrix  material  shrinks  upon  drying  of  the  composite. 
The  shrinkage  cracks  are  not  limited  in  the  matrix  but 
sometimes  extend  through  the  yarns.  The  micrographic 
images  in  Fig.  4  shows  these  cracks. 

d.  Interlaminar  separation:  This  type  of  defects  occurs  due  to 
a  low  level  of  compression  of  layers,  a  kinked  layer  or  a 
large  air/gas  pockets  that  may  occur  during  manufactur¬ 
ing.  They  show  in  micrographic  images  between  fabric 
perform  layers.  Fig.  5  shows  the  location  of  these  defects 
in  the  micrographic  images  and  in  side  and  planar  views. 

(3)  Architectural  defects  that  affects  the  preform  of  the  layer  as 

a  whole  such  as: 

a.  Layer  misalignment:  Misalignment  of  lamina  can  happen 
during  hand  layup  of  fabrics  due  to  error  in  placement 
of  preforms  prior  to  infiltration.  In  this  case,  the  cross 
section  of  the  fibers  in  the  micrographic  images  will 
not  appear  as  a  circle  but  rather  as  an  ellipse. 

b.  Fabric  wrinkling:  Where  the  fabric  may  suffer  wrinkling 
due  to  rough  handling. 

For  a  complete  understanding  of  the  3D  geometry  of  a  single  de¬ 
fect,  a  systematic  sectioning  of  the  composite  at  a  constant  spacing 
is  needed  to  trace  the  exact  morphology  of  the  defect.  But  it  also  is 
important  to  remember  that  many  defects  of  the  same  type  typi¬ 
cally  show  in  each  image  revealing  the  change  of  the  defect  geom¬ 
etry  along  its  path. 

5.  Analytical  model  of  elastic  properties 

A  model  was  constructed,  to  account  for  the  effect  of  the  above 
mentioned  defects  using  a  modified  stiffness  averaging  technique. 
Stiffness  averaging  was  selected  due  to  the  direct  application  of  the 
preform  geometry  in  its  formulation  and,  accordingly,  its  ability  to 
evaluate  the  effect  of  defects  at  different  locations  in  the  preform. 
The  model  takes  as  input  stiffness  of  various  constituents  of  the 
composite,  their  relative  volume  fraction  and  the  volume  fraction 
of  defects  in  different  directions.  A  set  of  analytical  relationships, 
detailed  in  [19],  use  these  parameters  to  evaluate  the  stiffness  ma¬ 
trix  for  each  local  yarn  direction.  These  stiffness  matrices  is  then 
transformed  to  the  material  global  directions  and  averaged  accord¬ 
ing  to  their  volumetric  contribution. 

To  ensure  the  overall  accuracy  of  this  approach,  elastic  proper¬ 
ties  calculated  using  this  model  for  a  hypothetical  non-defective 
composite  were  compared  to  those  calculated  using  the  hybrid  Fi¬ 
nite  Element  approach  (pcGlNA©)  reported  by  Gowayed  et  al. 
[9,10]  and  Gowayed  and  Barowski  [8]  and  the  difference  in  results 
for  both  models  was  within  5%. 
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Fig.  2.  Crossover  defects  in  a  micrographic  image  and  as  represented  using  pcGINA©.  The  dark  object  located  between  the  crimped  yarns  in  the  figure  on  the  left  represents 
this  type  of  defect. 


Fig.  3.  Matrix  voids  in  a  micrographic  image  and  as  represented  using  pcGlNA©.  The  dark  shapes  along  the  edges  of  the  yarn  in  the  figure  on  the  left  represent  this  type  of 
defect. 


The  stiffness  averaging  approach  divides  the  composite  into  a 
group  of  yarn  elements,  in  which  each  element  is  assumed  to  con¬ 
sist  of  thousands  of  coated  fibers  with  the  matrix  material  sur¬ 
rounding  each  fiber.  These  yarn  elements  have  lengths  and 
directions  that  satisfy  the  geometry  of  the  fabric  (i.e.,  volume  frac¬ 
tion  and  angle  of  inclination  for  each  yarn  segment).  Elastic  prop- 


Fig.  4.  shrinkage  cracks  in  oxide/oxide  composites.  Notice  the  almost  vertical 
cracks  enclosed  in  the  dotted  box  in  the  image. 


erties  of  the  composites  are  calculated  as  follows:  (i)  Calculate  the 
stiffness  properties  of  each  yarn  element  with  respect  to  its  local 
orthogonal  axis  system  (1,  2,  3)  in  which  axis  1  is  typically  parallel 
to  the  fiber/yarn  longitudinal  axis,  (ii)  rotate  the  local  stiffness 
matrices  in  space  to  evaluate  their  stiffness  with  reference  to  the 
composite  axis  system  (x,  y,  z),  in  which  in  this  work,  x  is  parallel 
to  the  warp  direction  (which  is  typically  oriented  along  the  length 
of  a  tensile  specimen)  and  y  is  parallel  to  the  fill  direction  (which  is 
typically  oriented  in  the  short  direction  of  a  tensile  specimen  in  the 
same  plane  as  the  warp  yarns),  and  (iii)  Average  the  stiffness  ma¬ 
trix  of  all  the  rotated  yarn  elements  with  respect  to  their  relative 
fiber  volume  fraction  in  the  composite.  Details  of  this  model  can 
be  found  in  the  work  by  Pastore  and  Gowayed  [19]. 

Although  each  category  of  defects  outlined  above  has  its  own 
directional  nature,  as  it  relates  to  the  fabric  architecture,  all  catego¬ 
ries  typically  show  large  dimensional  variability.  In  other  words, 
their  spatial  orientation  is  predictable,  but  their  spatial  dimension 
in  a  particular  direction  is  not.  Accordingly,  the  impact  of  defects 
on  elastic  properties  was  taken  as:  (i)  volumetric  by  the  reduction 
in  matrix  volume,  as  well  as,  (ii)  directional  by  including  the  effect 
of  the  projection  of  the  defect  along  each  direction. 

For  defects  with  a  clear  directional  nature  (i.e.,  dry  fibers,  cross¬ 
over  points,  shrinkage  cracks  and  interlaminar  separation),  two 
types  of  unidirectional  yarn  elements  are  identified:  a  non-defective 
and  a  defective  element.  The  elastic  properties  of  the  non-defective 
element  are  calculated  following  the  typical  procedure  of  the 
stiffness  averaging  technique,  while  those  for  the  defective  element 
are  calculated  based  on:  (i)  including  the  volume  of  defects  in  calcu¬ 
lating  the  volume  fraction  of  the  matrix  and  (ii)  using  the  vector  of 
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Fig.  5.  Interlaminar  separation  defects  in  a  micrographic  image  and  as  represented  using  pcGINA©.  The  dark  objects  located  above  the  crimped  yarn  in  the  figure  on  the  right 
represent  this  type  of  defect.  The  bottom  figure  show  possible  locations  of  this  defect  in  the  planar  view. 


the  area  projections  to  reduce  the  value  of  the  moduli.  Based  on  this 
approach,  the  equations  for  the  three  elastic  (£/d)  and  three  shear 
(Gyd)  moduli  (i,  j  =  1,  2,  3  and  i  for  defective  unidirectional  ele¬ 
ments  are  shown  in  Eq.  (1),  where  for  the  local  yarn  axis  system 
(1,  2,  3),  axis  1  is  in  the  longitudinal  direction  of  the  unidirectional 
yarn  element,  and  axes  2  and  3  are  orthogonal  to  1 : 

Eu=Ei{1-AVi) 

^  /1  AVj+AVj-. 


where  E/  and  Gy  {ij  =  1, 2, 3  and  i  ^j)  are  the  elastic  and  shear  mod¬ 
uli  of  the  non-defective  unidirectional  element,  and  Avi  (i  =  1,  2,  3) 
are  the  area  fraction  of  the  defects  on  a  plane  perpendicular  to  axis 

1  (i.e.,  area  of  defects  divided  by  the  total  area).  The  volume  of  de¬ 
fects  is  taken  into  consideration  by  reducing  the  volume  of  the 
matrix. 

The  vector  of  area  fraction  of  defects  (AVi)  will  change  from  one 
defect  type  to  another  based  on  the  geometric  characteristic  of 
each  defect.  Table  1  provides  schematic  drawings  for  the  defects 
and  the  equations  used  to  calculate  the  vector  of  area  fractions 
for  a  yarn  in  the  warp  direction.  For  the  local  axis  system  (1,2, 
3),  the  1  axis  is  in  along  the  fiber/yarn  longitudinal  axis,  and  axes 

2  and  3  are  orthogonal  to  1.  For  the  global  axis  system  (x,  y,  z) 
shown  in  Table  1,  the  x  axis  is  in  the  warp  direction,  the  y  axis  is 
in  the  fill  direction  and  z  axis  is  in  the  through-thickness  direction. 
Following  the  same  logic,  similar  equations  were  formulated  for 
the  fill  yarns  in  the  case  of  balanced  and  unbalanced  weaves. 

In  Table  1,  a,  h  =  long  and  short  diameters  of  the  yarn  elliptical 
cross-section,  respectively,  w  =  warp  and /=  fill/weft  directions. 
Pi  =  yarn  spacing  in  i  direction  (z  =  w,  /),  Ut  =  number  of  ends  in  i 
direction  per  unit  length  (z  =  w,/),  and  df  is  the  length  of  the  projec¬ 
tion  of  the  defect  in  the  fill  direction,  dw  is  the  length  of  the  projec¬ 
tion  of  the  defect  in  the  warp  direction,  dz  is  the  length  of  the 
projection  of  the  defect  in  the  through-thickness  direction  and  (*) 
is  used  as  a  multiplication  sign. 

For  the  dry  fibers  defect,  Eq.  (1 )  for  Fid  is  modified  to  include  the 
effect  of  the  stress-discontinuity  in  the  fibers  direction  utilizing  a 


shear-lag  model  as  follows,  while  the  rest  of  the  equations  for 
other  composite  moduli  remain  the  same: 


Eld 


EflSfVfTf  \ 
4t(L  -  dw)  J 


(l-Az^i) 


(2) 


where  Efi  is  the  fiber  longitudinal  modulus,  5/  the  fiber  strength,  Vf 
the  fiber  volume  fraction,  r/the  fiber  radius,  L  the  distance  between 
defects,  and  t  is  the  interfacial  shear  strength  between  the  fiber  and 
the  matrix. 

Matrix  voids  were  assumed  to  be  spherical  and  to  affect  the 
elastic  tensile  and  shear  moduli  of  the  matrix.  These  effects  were 
determined  utilizing  the  approach  developed  by  Weng  and  Tandon 
[21,22]  as  shown  in  Eq.  (3)  for  the  effect  of  inclusions  (r)  on  the 
bulk  {Kmv\  tensile  (Emv)  and  shear  (G^v)  moduli  of  the  homoge¬ 
neous  materials: 

Gm»  =  Cm  (l  + 

Emv  —  2Gf7it'(^  “1“  ^m)  fr>\ 


3km  +  4Gm 
_  ^{km  +  2Gm) 

“5(3/<„,+4G„,) 

For  a  matrix  with  inclusions  r: 


Vr(I(r-I(m) 

'  ^a(Kr-Km)+I<m 


L  _  W(Gr  —  Gm) 

’  “4"/^(Cr-G„)  +  G„ 

where  K,  G,  E  and  v  are  bulk,  shear  and  elastic  moduli  and  the  Pois¬ 
son’s  ratio  of  the  matrix,  respectively,  V  is  volume  fraction,  m  refers 
to  matrix,  r  refers  to  inclusion,  refers  to  voids,  a,  p,  and  bi  are  de¬ 
fined  constants.  The  inclusion  in  the  case  of  defects  is  assumed  to  be 
air. 
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Table  1 

Vector  of  area  fractions  for  different  defect  types  for  a  yarn  in  the  warp  direction. 
Defect  type  Schematic  of  defect 

Dry  fibers  *  z 


Crossover  points 


Shrinkage  cracks 


Interlaminar  separation 


r 

X 

dw 

Pf  =  1/  Hf 

Pf  =  1/  Hf 

Pf  =  1/  If 

X, 

z 

,a/2  -  b/2  ^a/2  -  b/j 

■  i 

^  Pf=l/nf  I 

2 

^  Pf=l/nf  1 

Pf  =  1/  rif 

x^ 

^  Z 

3= 

T3 

— 1 - 

^  1 
Pf  =  1/  n,  j 

^  Pf  =  1/  Hf 

Pf  =  1/ 

\  z 

X 

dw 

,  Pf  =  1/  nf 

1 

1  Pf=l/n, 

Pf  =  1/  Hf 

Area  fraction 


Av\ 

AV2 

AV3 


dz  df 
2bPw 
dz  dw 
2bpf 

dw  df 
PwPf 


Av\  = 
AV2  = 
AV3  = 


ba 

2bp„ 

^iPf-a+b) 

2bpf 

a(Pf-a+b) 

PwPf 


Av\  = 
AV2  = 
AV2,  = 


dz  df 
2bPw 
dz  dw 
2bp/ 
dw  df 
PwPf 


_ dw  df 

PwPf 


In  the  case  of  lamina  misalignment,  by  assuming  that  fibers 
have  a  circular  cross-section,  any  departure  from  such  geometric 
form  can  be  used  as  an  indicator  of  misalignment.  The  effect  of 
layer  misalignment  is  included  by  changing  the  orientation  angle 
in  the  rotation  of  the  stiffness  matrix  of  the  local  element  from 
the  local  axis  system  to  the  global  composite  axis  system  in  order 
to  match  the  actual  angle  of  the  layer. 

6.  Results 

Micrographic  images  obtained  from  different  samples  (10 
images  per  data  point)  were  thresholded  at  the  lower  infiection 
point  of  the  slope  of  the  grey-scale  curve.  Black  pixels  were  as¬ 
sumed  to  be  defects  and  were  manually  measured  and  categorized. 


Fig.  6.  Micrographic  (top)  and  processed  images  (bottom)  of  oxide/oxide 
composite. 


Instances  of  defects  were  categorized  per  defect  type.  For  each  de¬ 
fect  type  in  each  image,  different  instances  were  measured  and  re¬ 
sults  from  different  images  were  averaged.  Black  pixels  inside  the 
area  typically  occupied  by  the  yarn  and  running  parallel  to  the  fi¬ 
bers  were  assumed  to  be  fibers  pulled  out  during  polishing.  Figs.  6 
and  7  show  examples  of  image  processing  for  oxide/oxide  and  MI 
SiC/SiC  composites,  respectively.  These  thresholded  images  isolate 
the  defects  and  allow  the  quantification  of  their  area  or  length  frac¬ 
tions.  In  both  images,  the  large  ratio  of  the  length  of  the  interlam¬ 
inar  separation  defect  versus  the  total  length  of  the  image  can  be 
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Table  2 

Vector  of  area  fractions  and  volume  fraction  of  defects  (Vd)  for  a  yarn  in  the  warp  direction^  (in  percent). 


Intra-yarn 

Dry  fiber 

Inter-yarn 

Crossover  points 

Matrix  voids'^ 

Shrinkage  cracks 

Interlaminar  separation 

SiC/SiNC 

Av-i 

1.51 

37.8 

- 

- 

2.80 

AV2 

2.0 

9.9 

- 

- 

2.80 

AV3 

3.0 

29.9 

- 

- 

64.0 

Vd 

2.44 

4.51 

2.44 

- 

2.86 

Oxide/oxide 

Av-i 

0.60 

0.53 

- 

18.0 

1.55 

AV2 

0.72 

2.67 

- 

0.30 

1.55 

AV3 

1.2 

1.13 

- 

0.06 

31.9 

Vd 

0.79 

0.33 

0.79 

0.01 

1.57 

MI  SiC/SiC 

7\vi 

5.04 

3.15 

- 

- 

1.58 

AV2 

4.0 

1.24 

- 

- 

1.58 

AV3 

5.04 

3.12 

- 

- 

40.0 

Vd 

5.39 

0.4 

1.13 

- 

1.57 

^  For  the  balanced  fabrics  discussed  in  this  work,  vectors  of  area  fraction  in  the  warp  and  fill  directions  are  expected  to  be  similar. 
^  Matrix  voids  were  assumed  to  be  spherical. 


Table  3 

Properties  of  constituent  materials  at  room  temperature  as  used  in  the  model  (GPa). 


SiC/SiNC 

Oxide/oxide 

MI  SiC/SiC'" 

0.42^ 

0.46^ 

0.36^ 

Coated  fiber 

Efi 

184.7^" 

257.2'" 

332.37V 

Eft 

74.1 2^" 

257.2^ 

130.18V 

Ef 

31.20^ 

105.V 

60.05^ 

Vf 

0.16  V 

0.224^ 

0.178  V 

Matrix 

C-m 

139.70'^ 

212 

329.24^® 

Cm 

58.2V 

8.54^ 

139.27^ 

Vm 

0.20^ 

0.23® 

0.182® 

^  Data  from  manufacturer. 

^  Nano-indentation  [18]. 

Calculated  using  Rule  of  Mixtures. 

^  Calculated  using  elasticity  equations. 

^  Estimated. 

^  Uncoated  isotropic  fiber. 

®  Refs.  [7,11]  to  account  for  growth  of  necks  between  particles  via  surface  diffusion. 
Ref.  [9]. 


seen.  Table  2  lists  the  vector  of  area  fractions  and  volume  fraction 
of  each  defect  type  for  the  three  CMCs.  Shrinkage  cracks  were  lim¬ 
ited  to  oxide/oxide  and  no  layer  misalignment  was  observed  in 
these  images. 

Additional  data  needed  to  run  the  analytical  model  include 
information  on  the  preform  architecture,  detailed  earlier  in  the 
section  on  materials,  and  the  properties  of  the  coated  fiber  and 
the  matrix  as  listed  in  Table  3  at  room  temperature.  The  effect  of 
temperature  on  the  coated  fiber  for  the  SiC/SiNC  composite  is  listed 
in  Table  4.  Polymer  derived  SiC  and  SiNC  matrix  material  do  not  ex¬ 
hibit  a  major  change  in  their  elastic  properties  at  temperatures  be¬ 
low  1300  °C  [3],  accordingly,  the  matrix  properties  were  kept 
constant  at  different  temperatures. 


Table  4 

Properties  of  coated  SiC  fiber  in  SiC/SiNC  composite  at  elevated  temperatures  as  used 
in  modeP  (CPa). 


RT 

800  °C 

1100 °C 

1200 °C 

1300 °C 

Eft 

184.7 

171.7 

156.8 

145.5 

132.5 

Eft 

74.12 

72.4 

70.2 

68.4 

66.0 

Cf 

31.20 

30.5 

29.6 

28.8 

27.9 

Vf 

0.16 

0.16 

0.16 

0.16 

0.16 

^  Calculated  from  [2,16]  using  rule  of  mixtures  and  elasticity  equations.  Poisson’s 
ratio  was  assumed  to  remain  constant. 


In  Table  3,  %  are  the  tensile  modulus  of  the  fiber  in  the  lon¬ 
gitudinal  and  transverse  directions,  respectively,  G/  is  the  shear 
modulus  of  the  fiber  and  v/  is  the  Poisson’s  ratio  of  the  fiber,  Em 
is  the  tensile  modulus  of  the  matrix,  Gm  is  the  shear  modulus  of 
the  matrix,  is  the  Poisson’s  ratio  of  the  matrix  and  Vfis  the  fiber 
volume  fraction. 

Computer  codes  were  developed  utilizing  the  above  equations 
using  Matlab®  to  account  for  the  effect  of  each  defect  type  as  a  per¬ 
centage  of  depreciation  in  properties.  Utilizing  these  equations  and 
the  data  in  Tables  2-4,  the  effect  of  defects  on  the  in-plane  and  the 
through-thickness  moduli  as  well  as  the  in-plane  shear  modulus 
were  evaluated  for  the  plates  from  which  the  defect  data  was  ob¬ 
tained.  The  results  are  listed  in  Table  5,  with  at  least  5  samples 
per  data-point.  Layer  misalignment  was  not  observed  in  these 
composites  and  accordingly  their  results  are  not  reported.  Table  6 
shows  comparison  between  measured  data  and  calculated  proper¬ 
ties  with  the  effect  of  defects  mathematically  added.  The  effect  of 
temperature  on  the  in-plane  modulus  of  the  SiC/SiNC  composite 
is  shown  in  Fig.  8. 

7.  Analysis  and  discussion 

Utilizing  the  geometric  categorization  system  for  defects  out¬ 
lined  above,  work  done  and  reported  in  archived  literature  on  sim¬ 
ilar  composite  systems  can  be  classified  with  respect  to  the 
preform  architecture.  For  example,  the  work  done  by  Abul-Aziz 
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Table  5 

Analytical  results  showing  percent  reduction  in  elastic  property  due  to  defects. 


Intra-yarn  Inter-yarn 


Dry  fiber 

Crossover  points 

Matrix  voids 

Shrinkage  cracks 

Interlaminar  separation 

SiC/SiNC 

Ex 

2.80 

3.35 

3.94 

- 

3.51 

Fz 

1.36 

3.48 

3.64 

- 

63.36 

Oxide/oxide 

Ex 

1.84 

0.13 

0.44 

14.05 

1.36 

Ez 

0.13 

0.26 

1.44 

0.92 

31.93 

Gxy 

0.00 

0.07 

1.49 

8.80 

0.07 

MI  SiC/SiC 

Ex 

4.11 

0.51 

1.92 

- 

2.26 

Ez 

2.12 

0.00 

1.50 

- 

39.61 

Gxy 

1.20 

0.38 

1.42 

- 

0.65 

Where  E^,  are  the  tensile  modulus  of  the  composite  in  the  in-plane  and  through-thickness  directions,  respectively,  and  Gxy  is  the  in-plane  shear  modulus  of  the  composite. 

Table  6 

Experimental  data  and  analytical  results  (GPa). 

Experimental  data 

Analytical  results 

Difference  from  average  experimental  value  (%) 

SiC/SiNC 

Ex 

112.7  ±7.4 

112.8 

0.1 

Ez 

31.7  ±6.5 

29.3 

7.6 

Oxide/oxide 

Ex 

71.6  ±6.4 

68.5 

4.3 

Ez 

33.2  ±  8.8 

24.9 

25.0 

Gxy 

13.3  ±0.6 

13.3 

0.0 

MI  SiC/SiC 

Ex 

254.9  ±  7.4 

239.1 

6.2 

Ez 

130.0  ±3.4 

110.3 

15.2 

Gxy 

95.8  ±  0.6 

83.4 

12.9 

Fig.  8.  The  effect  of  temperature  on  the  in-plane  modulus  of  cross-ply  SiC/SiNC  8- 
layers  composite. 


et  al.  [1]  on  MI  SiC/SiC  composite  includes  matrix  voids  and  cross¬ 
over  points,  while  the  work  done  by  Peters  et  al.  [20]  on  CVI  SiC/SiC 
composite  only  considers  matrix  voids.  For  the  carbon/SiC  compos¬ 
ite  reported  by  Puglia  et  al.  [4,6],  the  A  classification  includes  dry 
fibers  and  matrix  voids,  the  B  and  C  classifications  are  similar  to 
shrinkage  cracks  and  the  D  classification  is  crossover  points.  No  re¬ 
cord  was  found  in  archived  literature  on  evaluating  the  effect  of 
interlaminar  separation,  dry  fibers  or  misaligned  layers. 

It  is  important  to  note  that  there  are  some  limitations  to  the  accu¬ 
racy  of  extracting  the  defect  data  from  the  micrographic  images.  For 
example,  interlaminar  separation  typically  occurs  at  the  edge  of 
yarns.  When  cut  along  the  length  of  the  yarn,  the  fibers  are  typically 
pulled  out  as  seen  in  Figs.  2  and  3,  thus  limiting  the  ability  to  identify 
the  exact  location  and  extent  of  the  interlaminar  separation.  Also,  in 
the  case  of  matrix  voids  and  dry  yarns,  the  resolution  of  the  images 
has  to  be  high  enough  to  allow  for  accurate  measurement  of  small 
amount  of  defects/voids.  The  issue  with  image  resolution  is  also 


crucial  when  it  comes  to  shrinkage  cracks  that  require  high  resolu¬ 
tion  as  well  as  high  magnification  in  order  to  capture  the  size  and 
distribution  of  the  cracks. 

Defect  types  can  have  different  distributions  in  different  com¬ 
posites  as  shown  in  Table  2.  Dry  fibers  were  highest  in  MI  SiC/SiC 
followed  by  SiC/SiNC  and  lowest  in  the  oxide/oxide.  Defects  at 
crossover  points  were  mainly  present  in  SiC/SiNC  composites.  Ma¬ 
trix  voids  were  highest  in  SiC/SiNC  and  much  lower  in  the  other 
two  composites.  The  effect  of  interlaminar  separation  was  more 
significant  when  considering  separation  area  rather  than  separa¬ 
tion  volume.  SiC/SiNC  had  the  highest  separation  area  followed 
by  the  MI  SiC/SiC  and  oxide/oxide  composites. 

Defects  were  prevalent  in  SiC/SiNC  composites  than  other  sys¬ 
tems.  This  is,  most  probably,  due  to  the  iterative  pyrolysis  steps 
used  during  manufacture  and  the  shrinkage  associated  with  each 
step.  Pyrolysis  of  preceramic  polymers  causes  them  to  remarkably 
shrink.  This  is  one  of  the  drawbacks  of  this  material  because  cracks, 
even  after  repeated  infiltrations,  remain  a  source  of  stress  concen¬ 
tration  in  the  matrix  as  well  as  a  source  for  reduction  in  thermal 
and  mechanical  properties.  The  shrinkage  is  due  to  two  reasons: 
(i)  Organic  and  inorganic  groups  oxidize  or  evaporate  as  well  as  a 
partial  reform  of  the  polymer  structure  due  to  crystallization, 
and  (ii)  At  a  very  high  temperature,  depending  on  the  polymer 
composition  and  the  pyrolysis  environment,  crystallization  starts 
and  causes  further  shrinkage,  causing  a  change  in  density  [13]. 

Defects  caused  reduction  in  the  elastic  properties  of  the  compos¬ 
ites  studied  in  this  work.  Interlaminar  separation  has  the  most  pro¬ 
nounced  effect  as  seen  in  Table  5  with  a  reduction  of  over  63%  in  the 
through-thickness  modulus  for  SiC/SiNC,  40%  for  the  MI  SiC/SiC  and 
around  32%  for  the  oxide/oxide.  Shrinkage  cracks  reduced  the  in¬ 
plane  modulus  of  the  oxide/oxide  by  about  14%  and  the  in-plane 
shear  by  8.8%  with  a  negligible  effect  on  the  through-thickness  mod¬ 
ulus  since  the  cracks  are  mostly  vertical  cracks  in  the  out-of-pane 
direction.  A  2.46%  of  matrix  voids  in  SiC/SiNC  composites  caused 
around  4%  reduction  in  the  in-plane  and  through-thickness  moduli. 
Dry  fibers  reduced  the  in-plane  modulus  by  4.1%  in  the  MI  SiC/SiC 
and  2.8%  in  the  SiC/SiNC  composites.  Also,  4.51%  crossover  point  de¬ 
fects  in  SiC/SiNC  composite  caused  a  reduction  in  the  in-plane  and 
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through-thickness  moduli  of  about  3.5%.  The  effect  of  other  defects 
was  equal-to  or  lower-than  2%  of  the  non-defective  elastic  proper¬ 
ties.  It  can  be  concluded  from  these  results  that  for  the  composite 
systems  investigated  in  this  work:  (i)  in-plane  modulus  is  mainly  af¬ 
fected  by  shrinkage  cracks  and  dry  fibers,  (ii)  through-thickness 
modulus  is  highly  affected  by  interlaminar  separation,  (iii)  in-plane 
shear  modulus  is  mostly  affected  by  shrinkage  cracks,  and  matrix 
voids,  and  (iv)  matrix  voids  and  defects  at  crossover  points  affected 
all  elastic  properties  by  an  amount  dependent  on  their  volume 
fraction. 

It  can  be  seen  from  Table  6  and  Fig.  8  that  model  predictions  are 
within  a  reasonable  range  from  experimental  data  for  most  proper¬ 
ties,  except  for  the  out-of-plane  properties  of  the  Ml  SiC/SiC  com¬ 
posite,  which  was  underestimated  by  around  15%  from  the 
experimental  range  of  the  data.  Interlaminar  defects  play  an 
important  dynamic  role  in  these  out-of-plane  properties.  The  re¬ 
gions  covered  by  these  defects  will  not  carry  tensile  forces  in  the 
out-of-plane  direction,  while  it  may  carry  compressive  stresses 
after  experiencing  a  limited  amount  of  deformation.  Since  the 
through-thickness  modulus  was  measured  by  applying  compres¬ 
sion  on  a  group  of  stacked  disks,  it  is  expected  that  its  value  in¬ 
creased  with  the  increase  in  the  stress  level  until  it  reaches  a 
final  plateau  which  can  be  as  high  as  the  value  of  a  defect-free 
composite.  A  detailed  discussion  of  the  through-thickness  com¬ 
pressive  modulus  of  Ml  SiC/SiC  can  be  found  in  [9]. 

Identifying,  categorizing  and  quantifying  defects  can  help  im¬ 
prove  the  composite  manufacturing  process.  For  example,  the 
shape  and  distribution  of  intra-  and  inter-yarn  defects  are  typically 
affected  by  manufacturing  pressure  and  temperature,  the  viscosity 
of  the  infiltrating  matrix  medium,  shrinkage  of  the  matrix  upon 
cooling  and  the  existence  of  rigid  particulates  as  part  of  the  matrix. 
A  study  of  the  effect  of  each  of  these  manufacturing  parameters  is 
needed  to  correlate  them  to  the  existence  and  distribution  of  a  spe¬ 
cific  defect  type  and  allow  for  their  reduction  and  eventual  elimi¬ 
nation.  For  example,  the  high  value  of  dry  fibers  in  Ml  SiC/SiC 
points  to  a  limited  penetration  of  the  CVl  SiC  into  the  yarns.  A  mod¬ 
ification  in  the  manufacturing  procedure  can  help  remedy  this  sit¬ 
uation.  Also,  shrinkage  cracks  are  only  evident  in  oxide/oxide 
because  no  additional  matrix  material  is  used  to  fill  up  these 
shrinkage  cracks.  In  the  case  of  SiC/SiNC  composite,  multiple  ma¬ 
trix  infiltration  steps  are  used  to  fill  up  voids  and  cracks  while  in 
Ml  SiC/SiC  Si  metal  infiltration  fills  up  the  voids. 

8.  Conclusions 

Defects  in  an  as-manufactured  oxide/oxide  and  two  non-oxide 
(SiC/SiNC  and  Ml  SiC/SiC)  ceramic  matrix  composites  were  catego¬ 
rized  and  their  volume  fraction  quantified  using  optical  imaging 
and  image  analysis.  It  was  observed  that  for  most  defect  types 
the  volume  fraction  of  defects  was  low  except  for  shrinkage  cracks 
in  oxide/oxide  composites,  crossover  defects  in  Ml  SiC/SiC  compos¬ 
ite  and  interlaminar  separation  for  all  composites. 

A  stiffness  averaging  approach  was  used  to  quantify  the  effect  of 
defects  on  elastic  properties  utilizing  the  value  of  their  volume 
fraction.  In-plane  tensile  and  shear  moduli  as  well  as  the 
through-thickness  compressive  modulus  for  the  composite  sys¬ 
tems  were  experimentally  evaluated  and  compared  to  results  of 
the  model.  It  was  observed  that  the  most  affected  property  was 


the  through-thickness  modulus  due  to  interlaminar  separation  de¬ 
fect  in  all  the  composite  systems.  In-plane  tensile  and  shear  moduli 
in  oxide/oxide  composite  were  also  reduced  by  shrinkage  cracks. 
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